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ABSTRACT 

This  report  summarizes  the  technical  work,  support, 
effort,  and  data  services  completed  during  the  Quarter 
January  through  March  1974.  The  report  is  organized 
into  Part  A  describing  Operations  and  Maintenance 
functions  (Tasks  1  through  4)  and  Part  B  describing 
Research  (Tasks  5  through  8).  In  Part  B,  we  describe 
only  work  represented  by  completed  technical  reports. 
Since  the  current  contract  replaced  formerly  separate 
contracts,  the  work  represented  by  one  of  the  technical 
reports  (SAAC  10)  was  accomplished  under  the  SAAC  con¬ 
tract.  Part  C  is  a  summary  of  the  energy  conservation 
measures  taken  at  the  SDAC  and  some  early  results  of 
those  measures.  Part  D  is  a  review  of  scientists  from 
foreign  nations  who  visited  the  SDAC  during  the  first 
quarter  of  1974. 
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A.  OPERATIONS  AND  MAINTENANCE  (Tasks  1  through  4) 


1 .  Operations 

During  the  quarter  LASA  data  was  recorded  93$, 

ALDA  97 o,  and  NORSAR  951  of  the  available  2160  hours. 

Most  of  the  outage  noted  for  LASA  was  caused  by  com¬ 
puter  failure  at  the  LDC.  Quality  checking  of  the 
NORSAR  transmission  line  produced  the  most  significant 
outage  for  this  link.  Transmission  error  rate  on  the 
TAL  for  this  report  period  was  1.0$. 

2 .  Programming 

a.  EP  and  DP 

During  this  quarter,  a  new  version  of  the 
DP  with  changes  to  reduce  low  rate  tape  error  was 
installed.  After  more  than  one  month  of  operation, 
the  error  rate  for  low  rate  tapes  was  reduced  by 
a  factor  of  20  to  1. 

Changes  were  made  to  the  EP  system  as  a 
result  of  the  elimination  of  the  E  and  F  rings  at 
LASA.  The  event  re-run  package  was  changed  to 
omit  use  of  E  and  F  rings.  Also,  the  EP  plot 
output  format  was  revised  to  account  for  the 
missing  subarrays  to  decrease  the  time  required 
for  EP  plot  preparation  and  plotting. 

b.  360/44 

The  OS  system  generation  was  started.  However, 
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some  difficult  problems  unique  to  the  360  Model 
44  have  been  found  and  we  are  working  on  them. 

Errors  affecting  the  machine  accounting 
reports  were  found  and  corrected  in  the  GRASP 
spooling  package. 

An  accounting  program  was  completed  during 
this  quarter.  It  is  expected  to  improve  uur 
360/44  utilization. 

We  have  installed  and  are  now  supporting 
PL1  as  an  application  language  available  on  the 
IBM  360/44  system. 

A  subroutine  was  completed  and  catalogued  to 
permit  full  use  of  the  commercial  instructions  on 
the  360/40's.  The  macros  "IF  THEN  ELSE"  and 
"DOWHILE"  were  catalogued  on  the  360/44.  They 
are  to  be  used  by  Network  Systems  to  support 
structured  programming  techniques. 

The  file  maintenance  program  was  modified 
to  accept  the  new  International  Seismic  Month 
data  tape  to  keep  our  capabilities  current  with 
the  development  of  this  data  base. 

c.  PDP-15 

The  PDP-15  required  engineering  changes  to 
shorten  memory  I/O  lines  and  the  I/O  bus  length. 
These  modifications  were  installed  by  DEC  factory 
engineers  and  supported  by  the  local  field  office. 
In  addition,  factory  representatives  from  BUCODE 
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installed  engineering  updates  to  the  9- track  tape 
units.  The  level  of  effort  rpplied  by  vendor 
support  amounted  to  about  40?o  of  prime  shift  time 

A  consultant  to  aid  our  staff  in  isolating 
operating  system  errors  was  obtained  during  this 
report  period,  tlis  assistance  was  used  during 
attempts  to  certify  the  hardware  and  system  soft¬ 
ware  modifications  made  throughout  this  quarter. 

SWAP  development  has  been  postponed  until 
the  system  hardware  and  software  modifications 
have  been  certified. 

The  effort  applied  to  the  Varian  software 
driver  permits  use  of  this  device  as  a  printer. 

d.  ARP A  Network 

A  program  has  been  written  for  the  360/91  to 
plot  seismic  data  via  network  protocol  on  our 
graphics  terminal.  The  purpose  of  this  program 
is  for  demonstration  and  to  help  us  understand 
network  graphics  conventions. 

We  have  completed  the  coding  and  testing  of 
our  network  software  to  transfer  seismic  data 
over  the  ARPA  Network.  Test  results  are  as 


follows : 

Stream  Mode  to  UCLA . . 14KB 

Block  Mode  to  UCLA - 1.4KB 

Stream  Mode  to  ILLIAC . . 2.0KB. 


The  stream  mode  transfer  rate  to  UCLA  is 
close  to  the  value  computed  by  theoretical  studies 
and  controlled  testing.  The  block  mode  transfer 
to  UCLA  produced  a  surprisingly  low  rate.  This 
low  rate  has  been  found  to  be  caused  by  the 
allocation  logic  in  UCLA’s  software.  We  have 
discussed  these  transfer  rates  with  their 
personnel.  Changes  to  their  system  must  be  made 
prior  to  further  testing. 

e.  Documentation 

The  following  documentation  was  completed 
and  delivered  during  the  quarter: 

36Q/Programs  PDP-15  Programs 

EVENT  CALCOMP  PLOT  Package 

M7APL0T  PLTPDP 

PLT360 

SWPXXX 

CLHAND 

Updates  to  ISRSPS  manuals  110S,  111S,  101, 
and  102  were  completed  and  delivered. 

3 .  Maintenance 

Line  voltage  monitor  equipment  for  the  PDP-15  was 
installed  for  one  month  during  this  report  period. 
Results  of  the  data  collected  during  this  test  are 
currently  under  investigation. 
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The  VLLA  time  code  reader  was  replaced  in  the 
analog  laboratory. 

The  magnifications  of  the  LASA  beam  devclocorder 
has  been  verified.  The  computed  and  measured  gains 
now  agree  at  76K. 

4 .  Data  Services 

a.  Analog  Processing 

Tape  compression,  A/0  conversion  and  special 
data  playouts  were  routinely  scheduled  after  we 
acquired  an  analog  technician  during  this  pciood. 
Currently,  30%  of  his  time  is  spent  on  A/O  con¬ 
versions,  and  30%  on  compression;  his  remaining 
effort  is  applied  toward  maintenance  and  satis¬ 
fying  special  data  requests. 

b.  Data  Recipients 

The  SDAC/LASA  weekly  event  summary  was 
mailed  to  thirty  individuals  or  institutions. 

The  following  individuals  and  institutions 
received  data  from  SDAC  during  the  quarter: 

MIT  Lincoln  Laboratory 

U.S.  Geological  Survey,  Menlo  Park 

Dr.  John  Lahr 
Dr.  William  Bakun 

Pennsylvania  State  University 

Dr.  Shelton  Alexander 
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U.  of  Texas  at  Dallas,  Dr.  Mark  Landisinan 
IJ.  of  Western  Ontario,  Mr.  Avadh  Ram 

Dr.  R.  1  .  Mereu 

HNSCO ,  Dr.  L .  Rage 
Institute  of  Geophysical  Sciences, 
Edinburgh,  Dr.  Graham  Neilson 
Texas  Instruments 

hlectronic  Signal  Processing,  Inc.  , 

Mr.  J.  Beardwood 

Epicenter  files  have  been  routinely  maintained 
as  the  data  became  available.  LASA  and  NORSAR 
bulletin  tapes  have  been  converted  to  a  format 
compatible  with  other  epicenter  files.  This  con¬ 
version  allows  a  search  to  be  made  of  the  various 
epicenter  files  for  common  events.  An  Earthquake 
Research  Laboratory  (formerly  NO A A ,  now  USGS) , 

LASA,  and  NORSAR  collation  file  is  routinely 
maintained.  As  of  the  end  of  the  quarter,  this 
collation  file  had  been  updated  through  30 
September  1973. 

c.  Tape  Libraries 

The  SDAC  seismic  data  library  contains  digital 
and  analog  magnetic  tape  and  film  recordings  from 
the  LASA,  NORSAR,  ALFA,  and  TFO  arrays.  It  also 
contains  analog  tape  and  film  from  the  seismic 
observatories  and  the  LRSM  sites  operated  early 
in  the  VLLA  program.  In  addition  to  the  primary 
data  described  above,  the  library  contains  digital 
tapes  either  converted  from  analog  data  or  refor¬ 
matted  from  other  digital  data  for  selected  events. 
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As  of  31  March,  there  were  approximately 
50,000  digital  tapes  and  56,500  analog  tapes  in 
our  library.  The  following  catagorics  of  digital 
tapes  are  permanently  retained  : 


LASA  I. vent  572 

NORSAR  Hvent  43 

Long  Period  (LASA,  NORSAR  4,058 

and  ALPA) 

Short  Period  LASA  3,508 

NORSAR  812 

Detection  Log  (International  30 

Seismic  Month) 

Lxtendcd  Long  Period  760 

UBO  Multiplexed  281 

LASA  Multiplexed  1,199 

TFO  Short  Period  649 

TFO  Long  Period  1,174 

A/D  and  D/A  Conversions  2,071 

Miscellaneous  1,137 


The  following  catagorics  of  digital  tapes 
are  recycled  or  used  in  operations  support: 

LASA  hi- rate  tapes  (LP  and  8,843 

SP  data) 

Detection  Log  14 

Individual  users  and  operations  2,800 

support 

The  analog  library  contains  the  following 
catagorics  of  tapes: 


Compressed  tapes  9,550 
Composite  tapes  332 
Tapes  as  recorded  18,000 
Scheduled  for  compression  8,925 


B.  RESEARCH  (Tasks  5  through  8) 


Six  reports  were  released  for  publication  during 
the  quarter.  One  of  these,  SAAC  10,  represents  work 
accomplished  under  a  former  contract;  the  remaining 
five  reports  concern  work  completed  under  the  SDAC 
contract.  A  summary  of  each  report  follows. 


1'  Detection  Threshold  of  the  LASA,  ALPA .  and  NORSAR 
Long- Period  Network  /SAAC  101 - 


Preliminary  estimates  of  the  network  thresholds 
for  the  detection  and  measurement  of  Rayleigh  waves 
were  obtained  in  this  study  using  long-period  data 
recorded  at  LASA,  NORSAR,  and  ALPA.  A  fast  frequency- 
wavenumber  detection  processor,  FKCOMB,  was  used  for 
analysis.  Events  used  in  this  study  were  obtained  from 
the  LASA/SAAC  Daily  Summary  and  the  NORSAR  Seismic 
Event  Summary  for  the  time  period  of  1  May  1971  to 
30  April  1972.  A  combined  ALPA-NORSAR-LASA  capability 
was  estimated  using  events  from  the  Kurile  Island 
Region,  and  a  combined  ALPA-NORSAR  capability  was  esti¬ 
mated  using  Eurasian  events.  For  the  ALPA-NORSAR-LASA 
network,  the  equivalent  body  wave  magnitudes  at  the 
90 »  level  of  detection  of  Rayleigh  waves  from  the 
Kuriles  are  as  follows:  ALPA,  mb  =  4.5;  NORSAR, 

mb  =  4,5»  LASA,  mb  =  4.5;  and  one  out  of  three  combined 
network,  m^  =  4.0.  From  cumulative  curves,  the  90% 
level  of  detection  of  Rayleigh  waves  from  the  Kuriles 
are  as  follows:  ALPA,  Mg  =  2.5;  NORSAR,  M  =  2.9;  and 
LASA,  Mg  =  3.0.  For  the  ALPA-NORSAR  network,  equivalent 
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body  wave  magnitudes  at  the  90%  level  of  uetection  ol 
Rayleigh  waves  from  Lurasia  are  as  follows:  ALFA, 
mb  =  4.7;  NORSAR,  mb  =  4.5;  ALPA- NORSAR  combined, 

=  4.4.  Prom  cumulative  curves,  the  9 0 '  level  of 
detection  of  Rayleigh  waves  from  Purasia  are  as  follows: 
ALPA,  M.  =  3.0;  and  NORSAR,  Mg  =  3.0.  Thresholds  were 
recomputed  for  the  respective  combined  network  wherein 
a  sample  was  not  deleted  if  a  Rayleigh  wave  was  masked 
at  one  array  and  was  not  detected  at  the  other  array. 

The  equivalent  mb  values  at  the  90%  threshold  were  as 
follows:  ALPA-NORSAR-LASA  for  Kuriles,  m^  =  4.2;  and 

ALPA-NORSAR  for  Eurasia,  mb  =  4.5.  The  thresholds 
calculated  for  ALPA  using  f-k  processing  agree  with 
those  established  using  other  methods  such  as  matched 
filtering  suggesting  strongly  that  thresholds  are 
geophysically  real  and  not  limited  by  processing  methods. 

2.  A  Comparison  of  the  Location  Refinement  Techniques 

In  the  SDAC/LASA  Event  Processor  fSl)AC-TR- 7$- 5) 

Two  techniques  for  refining  locations,  beampacking 
and  crosscorrelation,  which  are  programmed  into  the 
SDAC/LASA  automated  Event  Processor  (EP)  ,  are  shown  to 
produce  equivalent  array  beam  traces  for  LASA  short 
period  data.  The  crosscorrelation  method  requires  an 
operational  s ignal - to- noise  threshold  significantly 
above  the  SDAC/LASA  detection  threshold  of  lOdb.  In 
addition,  the  method  yields  unreliable  locations  for 
events  with  low  coherence,  such  as  events  with  low 
signal- to-noise  ratios,  mixed  events,  or  events  at 
short  epicentral  distances  from  the  array.  On  the 
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other  hand  the  beampacking  method  is  independent  of 
s ignal- to-noise  ratio.  Moreover,  beampacking  uses  15 
to  20  percent  less  computer  time  than  crosscorrelation 
and  obtains  locations  with  smaller  location  errors  when 
compared  to  Worldwide  Network  locations. 

This  report  also  shows  that  the  number  of  events 
reportable  on  the  SDAC/LASA  Daily  Summary  could  be 
increased  15  to  20  percent  by  lowering  the  12 P  acceptance 
threshold  from  14db  to  12db.  When  operating  the  HP 
with  a  s ignal- to-noise  acceptance  threshold  of  3  (lOdb) , 
beampacking  increased  the  number  of  verified  events  on 
the  Daily  Summary  by  approximately  30-40  percent  over 
the  number  produced  by  crosscorrelation. 

3.  LASA  Regional  Travel-Time  Corrections  and 

Associated  Nodes  (SDAC-TR- 73-6) 

This  report  describes  a  new  set  of  region  correc¬ 
tions  for  LASA  which  more  adequately  covers  the  seismi- 
cally  active  areas  of  the  earth  than  did  previous  sets. 
The  new  set  contains  183  calibration  nodes  versus  105 
nodes  on  the  set  used  in  the  SDAC/LASA  system  through¬ 
out  1972  and  the  first  half  of  1973.  Each  node  contains 
the  average  relative  travel-time  residuals  for  each  of 
the  21  subarrays  at  LASA  which  are  valid  over  an  area 
surrounding  the  location  of  the  node.  The  correction^ 
were  generated  from  more  than  1800  events  using  as  a 
primary  source  of  data  the  time  shifts  resulting  from 
the  cross c.orre lat ion  process  in  the  SDAC/LASA  Event 
Processor. 
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4. 


False  Alarm  Probabilities  for  Mixed  Lvents 

r^AC-Tk'-'75'-T) 


The  analysis  of  1,471  P-  and  PKP  codas  indicates 
that  the  probability  o f  an  unexplained  phase  occurring 
in  a  coda  of  an  event  as  recorded  at  a  single  station 
is  0.12  for  a  detection  threshold  on  the  order  of  3.Sdb 
(s ignal- to-coda  background).  The  average  coda 
length  is  roughly  G  minutes  (343  seconds)  for  the 
events  examined.  The  probability,  therefore,  that  the 
seismograms  at  four  stations  out  of  thirteen  will  ex¬ 
hibit  unexplained  phases  in  the  coda  of  a  given  earth¬ 
quake  is  0.045.  The  probability  that  unexplained 
phases  at  four  stations  will  yield  a  significant  loca¬ 
tion  solution  is  0.032.  Thus,  with  about  10,000  events 
occurring  each  year  of  magnitude  m^  >  4.0,  and  with  nil 
of  them  examined  for  unexplained  phases,  we  expect  15 
false  alarms  per  year.  The  probability  that  four  or 
more  stations  will  experience  a  false  alarm  is  essen¬ 
tially  the  same  as  the  probability  that  exactly  four 
stations  will  experience  a  false  alarm.  For  purposes 
of  on-site  inspection,  residual  travel-time  errors 

-3  <  c.  <5  seconds  imply  that  an  event  can  be  located 

1  5  2. 

anywhere  in  an  area  3.2  x  10  km  in  size. 

5 .  Seismic  Distance-Amplitude  Relations  for  Short- 
Period  P  ,  ff ,  PP  and  Compress ional  Core 
Phase's"  for '~A>  90°  (SbAC-tk-  73-9') - 

Measurements  of  log^(A/T)  reported  by  the  Inter¬ 
national  Seismic  Centre,  the  VLLA  Observatories,  and 
the  Long-Range  Seismic  Measurements  Program  are  used 
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in  this  study  to  define  an  ampl i tude- d is tance  curve 
tor  the  maximum  amplitude  in  the  first  few  seconds  of 
motion  for  the  distance  range  A>90°.  In  general  terms, 
the  corresponding  phases  are  P,  Pdifff  and  PKIKP. 

Some  information  is  also  obtained,  however,  for  the 
phases  PP  and  PKP^. 

0.  Regional  Attenuation  of  Short-Period  P  and  S 

Waves  in  the  United  States  ( SDAC-TR- 74 -  lT 

Regional  distribution  of  anelastic  attenuation 
beneath  the  United  States  was  investigated  using 
amplitudes  and  dominant  periods  of  short  period  P  and 
S  waves  originating  from  deep  focus  earthquakes  in 
South  America  and  the  Circumpacific  seismic  belt, 
recorded  at  LRSM  (Long  Range  Seismic  Measurement) 
stations.  The  observed  regional  distribution  pattern 
shows  high  attenuation  in  the  western  United  States, 
including  California,  and  a  less  pronounced  high 
attenuation  region  in  the  northeastern  United  States. 
This  distribution  pattern  is  similar  to  that  reported 
by  Solomon  and  Toksoz  (1970)  for  long  period  S  waves, 
but  differs  from  it  sufficiently  to  indicate  lateral 
variations  in  the  frequency  dependence  of  the  average 
crust-upper  mantle  attenuation  across  the  United  States. 
The  relative  changes  of  short  period  I’  and  S  wave 
amplitudes  can  be  sufficiently  explained  by  assuming 
a  complex  shear  modulus;  no  losses  in  compression  are 
indicated. 
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C.  ENERGY  CONSERVATION 

In  accordance  with  national  policy,  we  have  taken 
a  series  of  steps  to  conserve  energy  in  our  operations. 

We  have  reviewed  our  travel  requirements  to  determine 
which  trips  are  essential  and  to  assure  that  only  one 
person  travels  when  more  people  are  not  actually  re¬ 
quired  to  accomplish  a  given  mission.  In  the  parts 
of  our  facilities  where  a  special  environment  is  not 
required  for  computers,  we  have  taken  measures  to 
control  heating  and  air  conditioning.  Specifically, 
we  maintain  heating  thermostats  at  or  below  68°F  and 
cooling  thermostats  at  78°F;  outside  the  work  shifts 
we  turn  off  the  systems  entirely. 

We  have  made  two  studies  to  determine  the  effective¬ 
ness  of  our  efforts.  In  an  attempt  to  reduce  electrical 
consumption,  we  have  had  a  vigorous  campaign  to  reduce 
lighting,  the  results  of  which  are  shown  in  Figure  1. 

As  can  be  seen,  the  average  for  the  first  quarter  for 
1974  has  been  reduced  about  301  from  that  for  similar 
periods  in  1972  and  1973. 


After  a  campaign  to  reduce  the  use  of 
employees  coming  to  work,  we  have  produced 
following  results:  of  all  employees  in  the 
Laboratories  of  Geotech  381  participate  in 
99o  use  the  bus,  5%  use  motorcycles,  and  109o 
ride  bicycles  to  get  to  work. 
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D.  VISITING  SCIENTISTS 


During  the  quarter,  scientists  from  France, 

Russia,  and  Sweden  visited  the  Seismic  Data  Analysis 
Center,  The  results  of  these  visits  are  summarized 
below: 

Dr.  Guy  Kuster  and  Dr.  Bernard  Massinon  of  the 
Atomic  Lnergy  Commission  of  France  (CHA)  visited  the 
Seismic  Data  Anal/sis  Center  in  Alexandria,  Virginia 
from  22  through  24  January  1974.  We  discussed  with 
them  the  automatic  detection  algorithms  used  in  the 
SDAC/LASA  system  and  the  Seismic  Wave  Analysis  Program 
to  be  used  with  the  PDP- 15/50  system.  They  described 
their  organization,  the  seismic  network  in  France,  and 
their  plans  to  use  a  PDP-11/45  in  their  system. 
Unclassified  documents  were  exchanged,  and  they  expressed 
the  desire  to  formalize  an  exchange  of  data  between 
the  two  countries.  They  are  particularly  interested 
in  receiving  our  SDAC/LASA  Weekly  Event  Summary. 

Dr.  I.  P.  Passechnik  and  Dr.  V.  F.  Pisarenko  of 
the  Academy  of  Sciences,  U.S.S.R.,  visited  the  Seismic 
Data  Analysis  Center  (SDAC)  at  314  Montgomery  Street, 
Alexandria,  Virginia  from  04  to  08  February,  1974, 
inclusive,  to  exchange  technical  information  relative 
to  seismic  monitoring  of  underground  nuclear  explosions. 
Dr.  E.  Willis  and  Dr.  C.  F.  Romney  of  the  Defense 
Advanced  Research  Projects  Agency  were  official  hosts 
for  the  visit.  During  the  week  the  ARPA  representatives 
were  supported  in  formal  lectures  and  informal  discussions 
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by  the  scientific  staffs  of  Tcledyne  Geotech  and  Texas 
Instruments.  The  technical  discussions  concerned 
principally  the  detection,  location,  and  identification 
of  underground  nuclear  explosions  and  earthquakes  and 
related  subjects  such  as  time-series  analyses,  source 
studies,  and  seismic  instrumentation.  In  their  lectures 
the  Russian  scientists  described  seismology  in  the 
Soviet  Union,  earthquake  prediction,  spectral  analysis, 
the  nature  of  the  core-mantle  boundary  from  PcP/P 
amplitude  observations,  and  peaceful  uses  for  nuclear 
explosions. 

Dr.  Dahlman,  Mr.  Israelson,  and  Mrs.  Hornstrom  of 
the  Swedish  Research  Institute  of  National  Defense 
visited  the  SDAC  on  21  March  1974.  Dr.  Dahlman  is  a 
primary  advisor  on  test  ban  negotiations  and  Director 
of  the  Hagfors  Observatory;  he  is  assisted  by  Mr. 
Israelson.  Mrs.  Hornstrom  is  a  systems  programmer. 

The  primary  purpose  of  their  visit  was  to  see  the 
PDP-15  computer,  because  they  expect  delivery  in  May  of 
a  similar  system.  They  were  also  interested  in  our 
research  and  analysis  work.  During  the  morning  we 
discussed  the  role  of  the  SDAC  in  the  VELA  Program 
and  the  SDAC/LASA  system.  The  activities  of  the  morning 
terminated  with  a  tour  of  the  facility.  In  the  after¬ 
noon  Dr.  Dahlman  described  his  organization  and  some 
of  his  research,  and  Mrs.  Hornstrom  was  briefed  on  the 
PDP-15  system. 
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